Large— scale structure from superdiffusion in a driven dissipative system 
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A system far from equilibrium is characterized by unconventional many-body dynamical effects, 
which can lead to anomalous density fluctuations and mass transport. Interestingly, these structural 
and dynamic features often emerge simultaneously in driven dissipative systems. Here we seek an 
origin of their co-existence by numerical simulations of a two-dimensional driven granular gas. We 
reveal a causal link between superdiffusive transport and giant density fluctuations. The kinetic 
dissipation upon particle collisions depends on the relative velocity of colliding particles, and is 
responsible for the self-generated large-scale persistent directional motion of particles that underlies 
the link between structure and transport. This scenario is supported by a simple scaling argument. 

PACS numbers: 05.70. Ln, 45.70.-n, 66.10.cg 



Many materials of biological or industrial importance 
are driven far from thermodynamic equilibrium by an 
imposed energy flux, mediated by e.g. motor proteins in 
living cells, or boundary-induced flow of fluids or par- 
ticulate matter [H-Q. Given a mechanism for energy 
dissipation and sufficient relaxation time, such systems 
may reach a statistical steady state in which macro- 
scopic quantities remain constant. Unlike thermody- 
namic equilibrium, however, the underlying sequence 
of micro-states admits cyclical currents [l|, and conse- 
quently driven dissipative systems exhibit a richer vari- 
ety of structure formation than thermal systems 0, Q ■ 
The transport of matter within such systems can also 
be anomalous, such as when the mean-squared particle 
displacements grow faster than linearly in time. Exam- 
ples of such superdiffusive systems include dusty plas- 
mas Q|, intracellular transport turbulent fluids 
self-propelled particles Q and granular media [1, Q . 

Interparticlc interactions, even if just steric hindcr- 
ance, inevitably lead to some form of correlation in the 
motion of nearby particles. It is therefore expected that 
superdiffusion should have a measurable many-body con- 
sequence and a corresponding spatial signature. There 
is some suggestion of this in the aforementioned mate- 
rials: Superdiffusion and a diverging dynamic correla- 
tion length have been observed in both experiments on 
sheared frictional granular media |8| and simulations of 
frictionless particles [§] near the 'jamming' transition di- 
viding rigid from non-rigid packings; and hydrodynamic 
equations for self-propelled particles admit superdiffu- 
sion and long-range ordering Q. These examples sug- 
gest superdiffusion can be linked with long range spatial 
correlations, but it is not clear for which systems, if any, 
the link is causal. 

Here we investigate the link between anomalous mass 
transport and large scale structure formation in simula- 
tions of a model driven dissipative system, in which par- 
ticles are uniformly agitated on the single-particle level 
and dissipate energy through short-range pair interac- 



tions. This flux produces superdiffusive particle trans- 
port that can be attributed to the spontaneous forma- 
tion of convective currents [|| . We provide a simple the- 
ory predicting long-range structure formation as a conse- 
quence of the anomalous transport, and confirm the ex- 
pected small wavenumber divergence in the static struc- 
ture factor. As an intermediate calculation, we predict 
and observe giant number fluctuations in which varia- 
tions in the particle number A exceed the classical A 1 / 2 
expectation. This phenomenon was previously predicted 
for active nematics and observed in vibrated granular 
rods 11011 . and also claimed for spherical particle mono- 
layers , but failure to reach ergodicity cast doubt upon 
this latter case. We confirm both ergodicity and conver- 
gence with system size, providing unequivocal evidence 
that giant fluctuations can exist in non-equilibrium sys- 
tems of isotropic particles. 

Model. — Our system, inspired by vibrated granular 
monolayers Q, consists of radially symmetric parti- 
cles with short-range repulsive, dissipative interactions, 
driven by a homogeneous and isotropic, Langevin-like 
force noise. The particles are discs a with polydisperse 
diameters d a and equal mass density. Two discs a and /3 
interact with equal-and-opposite forces when their cen- 
ters are separated by a distance R a/3 < \{d a + d@). The 
interaction has a repulsive conservative component of 
magnitude / cons = ^[l-2R afj /(d a +d p )} acting along the 
line of centers, and a dissipative term f lss = rj (v Q — v' 3 ) 
which reduces the relative velocity. These interactions 
conserve momentum but dissipate (kinetic) energy. The 
driving term consists of a spatio-temporally uncorrelated 
fluctuating Gaussian force field £(r, t) that obeys white 
noise statistics, (£i( r ij ii)£j(r 2 , t 2 )) = T8 ij 5(vi-Y 2 )8(t 1 - 
t-2). This field does not conserve momentum locally, but 
it is imposed globally to ensure a fixed system centre of 
mass. Note that there is no frictional term with an im- 
plicit solvent or substrate. 

Discs with diameters uniformly distributed over the 
range [0.7(d), 1.3(d)] are randomly placed in an L x L 
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simulation cell with periodic boundaries, until the re- 
quired area fraction <f) = L~ 2 ^ a ir(d a /2) 2 has been 
achieved. Here we consider only the case <fi — 0-5, corre- 
sponding to intermediate densities far below the jamming 
transition for this polydispersity, </>j « 0.843 [12J. Rather 
than fix the force noise T, we instead choose T to give the 
required mean kinetic energy K{t) = ^m a [v a (t)] 2 in 
the steady state, with m a the mass of particle a. This 
r was interpolated from data generated by a series of 
calibration runs for each r\. 

To determine when a statistical steady state has been 
reached, we measure the two-time mean squared dis- 
placement (Ar 2 (t w ,t w + <)) = (|r(< w + t) - r(t w )\ 2 } for 
particle displacements between times t w and i w + 1. Sta- 
tionarity is assumed when (Ar 2 (£ w , t w + i)) ceases to 
vary with t w and time translational invariance has been 
achieved, i.e. {Ar 2 (t w , t w + t)) = (Ar 2 (i)). All quanti- 
ties are expressed in dimensionless forms after suitable 
scaling by the bare distance (d), time to — y/ (d)(m)//x, 
dissipation coefficient r/o — \/ (d) and energy Kq = 
(d)[i. For this study, all 9 combinations of K(oo)/K = 
5 x 1CT 3 , 5 x 1CT 4 and 5 x 1CT 5 and 77/770 « 0.04, 0.08 
and 0.24 were used. Snapshots are given in (l3j |. 

Results. — For all K and 77 investigated, the mean 
squared displacement exhibits 3 regimes (see Fig. [TJ: 
(i) A rapid initial growth for trajectories shorter than 
the particle size, (Ar 2 (t)) -C (d) 2 ; (ii) Superdiffusivc 
motion at intermediate times, (Ar 2 (t)) ~ t 1+a with 
< a < 1; (Hi) Normal diffusion (Ar 2 (i)) ~ t at late 
times. This latter regime moves to later times for larger 
system sizes L, and we infer it is a finite size effect 
and the true asymptotic behavior is the superdiffusive 
regime (ii). The crossover from (i) to (ii) is also evident 
in the distribution P(Ar 2 , t) of squared particle displace- 
ments Ar 2 over the lag time t. For each regime it is 
possible to collapse P(Ar 2 ,t) onto a single curve after 
scaling the axes by a power of t while preserving normal- 
ization, P{Ar 2 ,t) = t-»p(Ar 2 /t u ), as shown in Fig. □ 
For short times we consistently find v close to 2, suggest- 
ing we are approaching the expected ballistic regime with 
P{Ar 2 1 t) ~ p ball ([Ar/i] 2 ) and Ar/t the velocity. The 
collapsed curve can be fitted to a 'stretched' exponential 
p ba11 ~ exp{— A{Ar /t) 3 / 2 }, consistent with velocity dis- 
tributions in granular gases [3]. More relevant here are 
large times for which collapse is possible with v « 1.5, 
but now the master curve is, to good approximation, a 
Gaussian with no fat tail. Superdiffusion is therefore 
identified with the anomalous broadening of the whole 
distribution, i.e. the variance a (f) ~ t v ~ t 1+a , rather 
than large jumps by a small subpopulation of particles 
as observed in high density frictional packings [8(. 

Fitted values of the exponents a and v for different 
K and 77 are given in Table H These values scatter 
around v = 1 + a = 3/2, and while random errors for 
some points are not consistent with this value, we cannot 
rule out small systematic errors of « 0.1 — 0.2 resulting 
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FIG. 1: (Color online) (Top panel) Mean-squared displace- 
ment (Ar 2 (t)) for the given system sizes L, K/Ko = 5 x 10~ 4 
and 77/770 = 0.08. The dashed line has slope 1.5. (Lower left 
panel) Probability distribution function (PDF) of displace- 
ments for L = 150(d) at short times, scaled by t 1 ' 8 . The 
dashed line is oc exp{— yl(Ar) 1 ' 5 }. (Lower right panel) Same 
for larger times, scaled by t 1 ' 45 . The dashed line is Gaussian. 

from the limited scaling regimes used for fitting. We 
believe that the true asymptotic value is a = 1/2 and 
is a consequence of the global driving and momentum 
conservation. To see this, first recall that (Ar 2 (t)) is 
directly related to the velocity autocorrelation function 
R(t) = (v(0) • v(i)), where v(t) refers to the velocity of 
a tagged particle at time t, 

(Ar 2 (t))^(|r(t)-r(0)| 2 ) = 2 f ds (t - s)R(s) , (1) 

Jo 

a result that assumes only steady state @, [l5| ■ The su- 
perdiffusive case of interest here corresponds to a diver- 
gent lim t _ i . 00 J dsR(s), i.e. R(t) ~ t~ b with b < 1 (we 
ignore the marginal case 6=1). According to |T|), this 
corresponds to (Ar 2 (i)) ~ t 1+a with a = 1 — 6. Thus 
superdiffusion corresponds to an anomalous slow decay 
of the velocity autocorrelation function. 

For our 2-dimensional system, the value 6 = i is ex- 
pected based on simple scaling arguments as explained 
in [l6|. In brief, over a time scale r, the random force 
noise changes the momentum of each particle a by an 
amount Ap" ~ r 1 / 2 . These momentum fluctuations are 
smoothed by the short-range repulsive interaction and, 
since momentum is conserved, can only spread diffusively. 
Therefore on a time scale r the momentum fluctuation 
Ap a will be smoothed over a volume ~ t d I 2 in D di- 
mensions. Thus momentum fluctuations should scale as 
Ap Q divided by the number of particles over which the 
momentum has been shared, i.e. ~ r 1 / 2 /r £) / 2 ~ t -1 / 2 
when D — 2, so R(t) ~ t~ x l 2 and 6 = a = ^ as claimed. 

A direct corollary of the slow decay of R(t) is that par- 
ticles drift over arbitrarily long distances. Put precisely, 
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TABLE I: Exponents a, b, v, a and /3 from power law fits 
for (Ar 2 (t)} (MSD), collapse of P(Ar 2 ,i) (PDF), persistent 
directed motion (PDM; eqn. ©), giant number fluctuations 
SN/N 1/2 (GNF) and static structure factor S(q). L = 150(d) 
and the numbers in brackets gives the error in the last digit. 



Parameters 


MSD" PDF PDM GNF S(q) b 
~t 1+a ^t L - b ~ N a ~ q -P 


K/Ko = 5 x lCT d 




77/770 = 0.04 
77/770 = 0.08 
77/770 = 0.24 


0.4(1) 1.35(5) 0.4(1) 0.2(l) c 1.75(5) 
0.60(5) 1.35(5) 0.5(1) 0.3(l) c 1.75(5) 
0.60(5) 1.3(1) 0.5(1) 0.3(l) c 1.7(1) 


K/K = 5 x 10~ 4 




77/770 = 0.04 
77/770 = 0.08 
77/770 = 0.24 


0.45(5) 1.45(5) 0.5(1) 0.2(l) c 2.00(5) 
0.55(5) 1.45(5) 0.50(5) 0.3(l) c 1.9(1) 
0.65(5) 1.40(5) 0.6(1) 0.3(1) 1.7(1) 


K/Ko = 5 x 10- b 




77/770 = 0.04 
77/770 = 0.08 
77/770 = 0.24 d 


0.45(5) 1.50(5) 0.5(2) 0.2(l) c 2.0(1) 
0.50(5) 1.50(5) 0.5(1) 0.3(l) c 1.8(2) 
0.65(5) 1.50(5) 0.6(1) 0.3(1) 1.8(1) 



"Fits to (Ar 2 (t)) = [C/t 1+a + E/t]- 1 for large t. 
'Fits to S(q) = Fq~@ + G up to q corresponding to ss 1 — 3{d). 
c Not converged with system size (exponent increasing with L) 
d Used larger system size L = 200 (d). 



for b < 1 the integral of R(t) = (v(0) ■ v(t)) 
respect to the lag time t is unbounded, i.e. 



t~ b with 



dtR(t) = (v(0) • Ar(t)> ~ t 



l-b 



b < 1 



(2) 



with Ar(t) = r(t) — r(0) as before. Thus particles will 
tend to move arbitrarily far in the direction of their initial 
motion, in the statistical sense of (jSJ. We refer to this 
as persistent directed motion. We have directly measured 
the integral ^fy in our simulations and in all cases found 
power law growth, as demonstrated in Fig. [3] Fitted 
values of the exponent are given in Table U and are in all 
cases consistent with the predicted value 1 — b = 1/2 = a. 

Snapshots reveal that this single-particle quantity has 
a corresponding spatial signature. As evident in Fig. 02 
particle motion becomes correlated over larger distances 
when longer time intervals are considered. To quantify 
this effect, note that according to (|2J), particles will move 
on average a distance Sr ~ r 1_b on a time scale r. Using 
the same diffusive dispersion relation for momentum fluc- 
tuations as above, over the same time scale particle mo- 
tion will become correlated over a range £ ~ t 1 / 2 . Thus 
of the N ~ (P particles in a region of size £, a fraction 
SN/N ~ Sr/£ will leave or enter the region, generating 
number fluctuations of magnitude 



SN ~l 



N~ 



(3) 



As long as ([2} is not subject to a finite size cut-off, per- 
sistent directed motion will apply over arbitrarily large 
length and time scales. Thus for all region sizes £, there 
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FIG. 2: (Color online) (v(0) • Ar(t)) normalized by (d) and 
the mean velocity (v) for the system sizes L given in the key, 
K/K = 5 x 10~ 5 and 77/770 = 0.24. The dashed line has a 
slope of 0.6. 





(a) (b) 

FIG. 3: (a) Particle displacements over a time interval 
t/to ~ 120 for the same parameters as Fig. [2l where arrows 
denote direction of total displacement over this interval, and 
light (dark) discs correspond to large (small) displacements. 
(b) Same starting configuration for t/to ~ 1200. The mean 
displacement is ~ 1.7(d) in fa) and ~ 10.5(d) in (b). Larger, 
color figures available from [13| ]. 



will be a corresponding time scale r ~ £ 2 driving the fluc- 
tuations in (J3|). Note that although inspection of snap- 
shots suggests persistent motion breaks down primarily 
in low-density regions, such velocity-density couplings 
are not included in this simple theory. 

These number fluctuations are 'giant' when the expo- 
nent on the right-hand side of ([3]) is larger than 1/2. It is 
then straightforward to map fluctuations for large N to a 
divergence at small wavelengt h q in the static structure 
factor, S(q) - q -im-^)+D] rj7|_ An examp i e of both 

SN/N 1 / 2 and S(q) is given in Fig. [fl clearly demonstrat- 
ing giant fluctuations and a divergent S(q — > 0). The 
corresponding exponents for various parameters are given 
in Table U which should be compared to the prediction 
SN/N 1 / 2 - N 1 / 2 and S(q) - <T 2 for b = 1/2 and D = 2. 
The number fluctuations are susceptible to finite size ef- 
fects and we were only able to attain convergence for 



4 



the 2 points marked in the table. To ensure convergence 
with time, SN was measured following two procedures, 
one in which temporal averaging is performed before spa- 
tial averaging (PI), and a second in which the order of 
averaging is reversed (P2). These two measures agreed 
for all but the largest N as in Fig. @] indicating ergodic- 
ity ll| (we also checked that the intermediate scattering 
function decayed by at least an order of magnitude over 
the same interval). Where available the fitted exponents 
are not inconsistent with the prediction, and we note the 
slope of the curves monotonically increases with system 
size in all cases, so giant fluctuations will become more 
pronounced, not less, for infinite systems. 

For S(q) the picture is clearer. Given a similar mag- 
nitude of systematic error induced by the narrow fit- 
ting regime as before, the S(q) exponents in Table U are 
within reasonable distance of the prediction. A diverging 
S(q) had also been observed in randomly driven inelas- 
tic hard-sphere systems, with the exponent 2 predicted 
by granular hydrodynamics [13, EH , suggesting its exis- 
tence may not depend on the details of the dissipation 
mechanism, i.e. whether it is scalar or vector. Super- 
ballistic mass transport was also claimed in hard-sphere 
simulations [Hj, but no causal link with S(q) suggested. 
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FIG. 4: (Color online) Number fluctuations 8N/N 1/2 ver- 
sus TV for different system sizes L, and for the two different 
procedures PI and P2 corresponding to time-averaging first 
(PI, solid lines) and spatial-averaging first (P2, dashed lines). 
K/Ko = 5xl0 -5 , 77 /770 = 0.24, and the dashed line has a slope 
of 0.31. (Inset) Static structure factor S(q) for the same sys- 
tems. The solid line corresponds to S(q) oc g~ 1,8 + B. 

We have thus demonstrated the co-existence of su- 
perdiffusion and large-scale structure in a driven dissi- 
pative system, and provided a simple theory implying a 
causal link between the two which agrees with the avail- 
able numerical data. We do not claim this relationship is 
general, and cite the contrary examples of turbulence in 
incompressible fluids, which can be superdiffusive with- 



out density fluctuations [6|, and isochoric critical fluids, 
which have a divergent S(q) but no superdiffusion [l9j . 
It is possible that superdiffusion inevitably leads to some 
form of long-range static or dynamic correlations that 
need not take the form of density fluctuations as in our 
model, but to confirm this would require more careful 
inspection of candidate systems. Additionally, our scal- 
ing theory does not include K or rj; empirically we find 
a non-trivial dependency of the various prefactors on 
these quantities which cannot be decomposed into sepa- 
rate power laws for each variable. An enhanced theory 
predicting the full scaling laws would be desirable, as 
would any attempt to deepen our understanding of this 
potentially far-reaching non-equilibrium phenomenon. 



[1] 

[2] 
[3] 



[4] 
[5] 

[6] 

[7] 

[8] 

[9] 
[10] 

[11] 



[12] 
[13] 



[14] 
[15] 
[16] 
[17] 
[18] 
[19] 



See e.g. D. Mukamel in Soft and Fragile Matter, eds. M. 

E. Cates and M. R. Evans (IOP, Bristol, UK, 2000). 
I. S. Aranson et al, Rev. Mod. Phys. 78, 641 (2006). 

T. Surrey et al., Science 292, 1167 (2001); T. B. Liver- 
pool et al, Phys. Rev. Lett. 90, 138102 (2003); K. Kruse 
et al, Phys. Rev. Lett. 92, 078101 (2004). 
S. Ratynskaia et al, Phys. Rev. Lett. 96, 105010 (2006); 

B. Liu and J. Goree, Phys. Rev. Lett. 100, 055003 (2008). 
P. Bursac et al., Nature Materials 4, 557 (2005); E. H. 
Zhou et al, Proc. Nat. Acad. Sci. 106, 10632 (2009); L. 
Bruno et al, Phys. Rev. E 80, 011912 (2009). 

A. J. Majda and P. R. Kramer, Phys. Rep. 314, 237 
(1999). 

Y. Tu, J. Toner and M. Ulm, Phys. Rev. Lett. 80, 4819 
(1998); H. Chate, F. Ginelli, G. Gregoire and F. Ray- 
naud, Phys. Rev. E 77, 046113 (2008); A. Baskaran and 
M. C. Marchetti, Phys. Rev. Lett. 101, 268101 (2008). 

F. Lechenault, O. Dauchot, G. Birolu and J. P. 
Bouchaud, Euro. Phys. Lett. 83, 46003 (2008); F. 
Lechenault, R. Candelier, O. Dauchot, J.-P. Bouchaud 
and G. Biroli, pre-print cond-mat/1001.1765. 

C. Heussinger, L. Berthier and J.-L. Barrat, pre-print 
cond-mat/ 1001.0914. 

S. Ramaswamy, R. A. Simha and J. Toner, Europhys. 
Lett. 62, 196 (2003); V. Narayan, S. Ramaswamy and N. 
Menon, Science 317, 105 (2007). 

I. S. Aranson, A. Snezhko, J. S. Olafsen and J. S. Urbach, 
Science 320, 612c (2008); V. Narayan, S. Ramaswamy 
and N. Menon, Science 320, 612d (2008). 

D. A. Head, Phys. Rev. Lett. 102, 138001 (2009). 

See EPAPS Document No. [] for snapshots and a larger, 
color version of Fig. [3] For more information on EPAPS, 
http:/ /www. aip.org/pubservs/epaps. html 



see 
T. 



CTvan Noije et al, Gran. Matt. 1, 57 (1998). 
G. I. Taylor, Proc. Lond. Math. Soc. 20, 196 (1922). 
A. Fiege et al, Phys. Rev. Lett. 102, 098001 (2009). 
T. P. C. van Noije et al, Phys. Rev. E 59, 4326 (1999). 
G. Peng and T. Ohta, Phys. Rev. E 58, 4737 (1998). 
A. Onuki, Phase Transition Dynamics (CUP, Cambridge, 
2002); J. K. Tison et al, J. Chem. Phys. 54, 1526 (1971). 



